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For slowly spinning matter the rate of energy loss via radiation of gravitational waves is estimated in 
General Relativity (GR) within a generally covariant superenergy approach. This estimation differs 
from Einstein's Quadrupole Formula (EQF) by a suppression factor (n <C 1). For a symmetric 
two-body- like distribution of scalar matter II is estimated to be <C (v / c) 2 (r / R) 2 , where v is orbital 
velocity of the bodies, c - velocity of light, r - radius of each body, and R - the inter-body distance. 
This contradiction with EQF is briefly discussed. 



Until the direct detection of the gravitational waves 
(GW) produced by binary neutron stars or other astro- 
physical and cosmological objects (which is expected as 
highly probable in the nearest decade) , the measurements 
of orbital damping of the binary pulsar PSR1913+16 
[1,2], in excellent agreement with Einstein's Quadrupole 
Formula (EQF) [3], are usually considered the ultimate 
indirect test of the energy radiated in General Relativity 
(GR) [4] by any moving matter. 

The generation of gravitational radiation in GR is a 
long-standing problem that dates back to the first years 
following the publications of GR [3,5]. There followed 
a lengthy debate about whether gravitational waves are 
real or an artifact of general coordinate invariancc. The 
fact that GW are real has been confirmed by coordinate 
free theorems [6-8] and by short-wave analysis [9] . 

Quantitative estimations of energy radiation via GW 
are based on (or agree with) the well known Einstein's 
Quadrupole Formula which is derived for pseudoenergy* 
radiation rate (pE) . Then the generally covariant energy 
radiation rate (E) is estimated by the noncovariant pE in 
the transverse-traceless (TT) gauge. The leading order 
of the post-Newtonian (PN) expansion yields 



*In GR the gravitational field (Levi-Civita connection co- 
efficients V kl ) does not possess any energy-momentum ten- 
sor but, as a consequence of the Einstein Equivalence Princi- 
ple (EEP), it only possesses the so-called "energy-momentum 
pseudotensors" (t M "). The gravitational energy-momentum 
(and gravitational angular momentum) pseudotensors, as be- 
ing functions of T kl (the total gravitational strengths) de- 
scribe the energy-momentum of the total gravitational field, 
which is a combination of the real gravitational field (for 
which the Riemann curvature tensor Rikim ^ 0) and the in- 
ertial forces field (for which Rikim = 0). The inertial forces 
field is generated by the coordinates used. This is also a con- 
sequence of the EEP. 



where G is Newton's gravitational constant^, and 
D ik = [ T 00 (3x t x k - 8 tk r 2 ) 



(1) 



(2) 



is the time-dependent, untraced mass quadrupole tensor 
with mass density T 00 given by the energy-momentum 
tensor of matter T^ v [10]. 

In spite of the generally noncovariant relation E = pE, 
after the discovery of binary pulsar PSR1913+16 [1], for- 
mula (1) with next-to-leading order PN contributions 
(see [11-13], in particular) is widely accepted for GR es- 
timations of GW radiation, in part because of theoretical 
works designed to shore up its foundations [12,13,15-20] , 
but mostly because of its excellent agreement with binary 
pulsar data [2,21]. 

But an empirical agreement does not imply conceptual 
adequacy; the covariance in theoretical derivations moti- 
vated by many experiments may be more important for 
the theoretical understanding of the nature of gravitation 
than a particular agreement with data. 

Gauge invariance is an important property of the wave 
equation and energy-momentum tensor for gravitational 
waves in GR. The usual gauge transformations used for 
gravitational waves are related to coordinate transfor- 
mations and gauge invariance implies invariance under 
these transformations. In general the wave equation for 
weak gravitational waves is approximately gauge invari- 
ant only for high frequency waves [9]. The effective 
energy-momentum tensor is gauge invariant to leading 
order only if it is averaged over a region of spacetime 
whose scale is large compared to the wavelengths of the 
waves. In other cases, such as low frequency gravitational 



t Throughout this paper the units, in which the light velocity 
i = 1 and reduced Planck constant h = 1, are adopted. 
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waves where the averaging is done over time, either the 
wave equation, the energy- momentum tensor, or both, 
are in general not gauge invariant [22]. 

Another alternative is to define gauge invariant effec- 
tive energy-momentum tensors for gravitational waves 
and other gravitational perturbations in almost all sit- 
uations of interest [22] if only those gauge transforma- 
tions are used which change the perturbed geometry, 
but leave the background geometry alone. This point 
of view allows considering the general covariance for all 
approaches to energy radiation based on the so called 
gauge invariant variables in GR [9,14,23-26]. Attempts 
to quantify the backreaction effects of GW in GR by 
defining a gauge invariant effective energy-momentum 
tensor for the waves were first done in [14] and later in 
[9] for high frequency gravitational waves in a vacuum. 
Another effective energy-momentum tensor for high fre- 
quency waves, both in a vacuum and in spacetimes con- 
taining classical matter, was defined in [23]. An extension 
of [9] in order to include lower frequency waves and space- 
times containing classical matter has been done in [25]. 
It should be noted that in each case mentioned above a 
procedure violating general gauge invariance was used: 
either some sort of averaging procedure, or a procedure 
of reduction of the gauge group that gives similar results 
to an averaging procedure [22] . 

In many known calculations with GW, the background 
spacetime manifold is considered non dynamical and, 
moreover, flat. This means a partial gauge fixing in order 
to work only with systems of coordinates such that the 
background metric takes the Minkowski form. What is 
left from the general gauge group after this gauge fixing 
is just the Poincare group. Thus, the reduction of the 
gauge group means the transition from a general invari- 
ant theory of gravitation to a Poincare invariant one [27] 
which is a theory of tensor gravitational field. 

This reduction of the gauge group is a possible rea- 
son why some formally independent calculations of GW 
radiation in the TT gauge ( [28-30] and [26]*) esti- 
mate the rate of energy loss in agreement with each 
other and with EQF: the energy-momentum pseudoten- 
sor transforms like the tensor under transformations of 
the Poincare group, as do the gauge invariant variables 
introduced with the partial fixing of the general gauge 
group. Namely, because of the partial fixing of the gauge 
group, all these estimations correspond to the tensor field 
gravity in Minkowski space-time and therefore they are 
close to each other. Unfortunately, all these estimations 
do not correspond exactly to GR. 

A classic point of view [31-33] is that the Riemann 
curvature tensor (R^ap) has to play the main role in 



'In these numerical calculations the gauge conditions are 
approximately transverse and traceless in the wave zone. 



the definition of gravitational radiation for the exact so- 
lutions to GR: only the covariant expressions dependent 
on R^uap may be used to get real information on grav- 
itational energy-momentum and angular momentum in 
arbitrary admissible coordinates. 

The canonical superenergy [34] and angular super- 
momentum [35] tensors are exactly this type of quan- 
tities. They extract covariant information (hidden in the 
pseudotensor (t^) [36]) about the real gravitational field, 
and can be expressed through d a pt^ v <~ Rz, va p [37]. 

Calculations of the superenergy radiation rate (sE) can 
be made by using conservation of the total superenergy 
introduced for a massive scalar matter field (4>) with met- 
ric gravitation [38]. 

As a contribution to debates in [17], [39], [40], [41,42], 
[43], [44], this paper is an attempt to consider the prob- 
lem of energy radiation by scalar matter in GR within 
a generally covariant superenergy approach [38] with the 
conservation of total superenergy and a covariant rela- 
tion between the energy-momentum and superencrgy- 
supermomentum tensors. 

Let's start from an approximate two-body-like solution 
to the Einstein-Klein-Gordon equation^ 

Rp* = 8ttG + \m 2 (t> 2 g^ , (3) 

where m is the scalar field mass. 

For semi-quantitative estimations let's approximate so- 
lutions to (3) by an extended, two-body-likc distribution 
of the scalar matter field <p(t, x, y, z) chosen in a form of 
fuzzy surface ellipsoid (with major semi-axis A, B, C) 
which is slowly spinning at angular velocity u>. In or- 
der to study the sensitivity of thus obtained results to 
deviation of <j) from the exact solution to (3), let's con- 
sider two different coordinate dependences for (f>(t, x, y, z) 
chosen in exponential <j>E(t, x, y, z) or Gauss <j>G{t, x, y, z) 
forms**. Neglecting gravitational mass shift, let's nor- 
malize both distributions to the rest mass of the matter 



§ Following [37,40], the convention is for Greek indices to 
run over four time-space values 0, 1, 2, 3, while Latin indices 
run over three spatial values 1, 2, 3; commas and d denote 
partial derivatives with respect to a chosen coordinate system, 
while semicolons and V denote covariant derivatives; repeated 
indices are summed over. 

**In the x — y rotation plane 4>e ~ 4>i, 4>G ~ 02, where 

^ j£sin(ag) -ycos(a;t)| ^ n ^— ^ j (4) 
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(M) at w = 0. For thus denned <pE and <f>a the energy 
and superenergy densities of the matter at A >• B <~ C 
are two-body-like (Fig.l). Although, ^ and 0g neglect 
radial contractions of each body under the gravitational 
attraction, any radial movement of matter does not ra- 
diate the tensor metric waves and is beyond the subject 
of this paper. 




FIG. 1. Two-body-like energy To (left) and superenergy 
Wo (right) densities projected on the rotation plane for (4e at 
u = 0, and A : B : C = 10 : 1 : 1. 

In the TT gauge small deviations (h al3 ) of metric ten- 
sor g^vix) from its Minkowski form {rj a P) yield [10] 
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(5) 



Instead of the noncovariant calculations with the 
energy-momentum pseudotensor and gauge invariant 
variables, let's calculate the evolution of the scalar mat- 
ter (S) [38] and metric gravitation (B) [34] terms of the 
conserving' 1 " 1 ' covariant superenergy tensor** (W a ^xn) [38] 



(7) 



In weak field approximation the calculations can be 
significantly simplified §§ , because the generally covariant 
superenergy density (W) and flux (Pj) can be chosen as 
W = VFoooo, and Pi = VF i0 oo- 

In order to compare rates of the superenergy and pseu- 
doenergy losses, let's substitute <j> E or (j> G into (5). This 



**A triple contraction of W a /3\^ with any three Killing vec- 
tors ( or three copies of the same if only one Killing 
vector exists) satisfy the differential conservation law [38] 



W a p X ^£%] 



0. 



(6) 



"Because the exact conservation of W (6) takes place only 
for the metric tensor and matter field satisfying (3) , all calcu- 
lations should not be too sensitive to a reasonable variation 
of the approximate solutions mentioned above. 

^For the goal of this paper the most reasonable choice of 
Killing vectors is the 4- velocity u of an observer at rest relative 
to the radiating matter: u = (1, 0, 0, 0). 



yields simple analytical expressions for corresponding 
fluxes. At large distances from the source of GW these 
fluxes are well known as the Bel-Robinson supermomen- 



tum (P -> P 



E.G 



[34] and Landau-Lifshitz pseudomo- 



mentum (Pj^) [10]. 

In the two-body-like case {A ^> B = C) calculations 
of the z-components of the fluxes generated by matter 
spinning in the x-y plane yield 



r ,GM 2 a4 „ 
R = A—^A 4 ^ 



r „GM 2 t4 o 
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(8) 



(9) 



Expression (8) differs from (9) by the factor 4w 2 , and the 
superenergy conservation law (6) results in a covariant 
quadrupole formula which differs from (1) by 4w 2 : 



W ~ 4w' 



G /■■■ik\ 

45> ) 



(10) 



Finally, E and E can be covariantly related to W and 
W. Substitution of (t> E or tfic into T 00 and VF 0000 re- 
lates the total energy (Ee,g) an d superenergy (We,g)- 
Neglecting the small high order terms (u 4 A 4 <~ v 4 <C 
J 2 A 2 ~ v 2 ) in the two-body-like case (A > B = C) 
yields the following covariant relations: 



E f 



2m 2 



We + \m, E g = + ^M. (11) 



The reasonably small difference between Ee and Eq 
roughly approximates possible deviations from the exact 
result. Because these deviations are not significant for 
semi-quantitative estimations, only the exponential case 
is considered below: 



E 



2m 2 



W. 



(12) 



Expressions (8), (9), and (12) together yield the gen- 
erally covariant quadrupole formula 



E 



G_ 

45 



n(zf) 



which differs from (1) by a suppression factor 



n«2— , 



(13) 



(14) 



where m (the scalar field mass) is yet an arbitrary pa- 
rameter. 

To be consistent with the classical solutions to the 
problem of GW radiation, m should be large enough: 
the scalar particle Compton wavelength (1/m) should be 
<C A,B,C, and correspondingly m 3> m , where the 
energy of quantum fluctuations mo is estimated by the 
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inverse values of A, B, C . Therefore, in accordance with 
(14) the upper limit of the energy loss rate can be es- 
timated roughly at m — to . The following particular 
choice too = V^4~ 2 + B~ 2 + C~ 2 significantly simplifies 
all analytical expressions. 

Taking into account that to » m , II is estimated to be 
extremely small in the two-body-like case: II <C (lot) 2 — 
v 2 (r/R) 2 <C 1, where r and R are characteristic scales 
(radial (R w A) and transverse (r « B w C)) of the 
matter distribution projected on the rotation plane, and 
v (= uR <C 1) is the orbital velocity. 

The kinematic factor v 2 in II results in that the quad- 
rupole gravitational radiation occurs at the order of v 7 
beyond Newtonian gravity instead of w 5 in the other ap- 
proaches with the partial violations of general covariance. 
The geometric factor (r/R) 2 suppresses the energy radi- 
ation to zero for point-like masses (r — > 0) in accordance 
with the equations for geodesies. 

Tipical values of II are estimated to be II <C 10~ 2 for 
the final stage of binary neutron stars (v <~ 0.3, r/R <~ 
0.5), and II <C 10~ 16 for the current stage of the binary 
pulsar PSR1913+16. These estimations for II strongly 
contradict to II = 1 in (1), and the binary pulsar data. 

A possible reason for the dramatic difference between 
the covariant (13) and noncovariant (1) estimations of 
the energy loss rates is the covariance of the superen- 
ergy approach which allows extracting only the covariant 
contributions from the energy-momentum pscudotensor 
. This pscudotensor contains contributions of iner- 
tial forces (with R a /3^ — 0) which may dominate in 
in some particular gauges. Nevertheless, these inertial 
forces do not contribute to the covariant energy flux of 
the real gravitational fields with R a /3^ ^ 0. The strong 
contradiction between the covariant and noncovariant 
calculations may mean far-reaching consequences for the 
structure of any metric theory of gravitation which claims 
to be realistic. 
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